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We present Monte Carlo simulations of surface induced disordering in a model of a binary alloy 
on a bcc lattice which undergoes a first order bulk transition from the ordered D03 phase to the 
disordered A2 phase. The data are analyzed in terms of an effective interface Hamiltonian for a 
system with several order parameters in the framework of the linear renormalization approach due 
to Brezin, Halperin and Leibler. We show that the model provides a good description of the system 
in the vicinity of the interface. In particular, we recover the logarithmic divergence of the thickness 
of the disordered layer as the bulk transition is approached, we calculate the critical behavior of the 
maxima of the layer susceptibilities, and demonstrate that it is in reasonable agreement with the 
simulation data. Directly at the (110) surface, the theory predicts that all order parameters vanish 
continuously at the surface with a nonuniversal, but common critical exponent (5\ . However, we find 
different exponents /3i for the order parameter (tp2,ip3) of the D03 phase and the order parameter 
ipi of the B2 phase. Using the effective interface model, we derive the finite size scaling function for 
the surface order parameter and show that the theory accounts well for the finite size behavior of 
{ip2, 1P3), but not for that of ipi. The situation is even more complicated in the neighborhood of the 
(100) surface, due to the presence of an ordering field which couples to ipi. 



I. INTRODUCTION 

First order phase transitions in the bulk of systems 
can drive a variety of interesting wetting phenomena at 
their surfaces and interfaces. They have attracted much 
attention over many years jjj, and are still very actively 
investigated Prominent examples are the wetting of 
a liquid on a solid substrate at liquid-vapour coexistence, 
or the wetting of one component of a binary fluid below 
the demixing temperature on the walls of a container. 
These systems are representatives of a generic situation, 
which has been studied in particular detail: Three phases 
coexist, substrate, liquid and vapour. The substrate acts 
as inert "spectator" which basically provides the "bound- 
ary conditions" for the liquid-vapour system. The liquid- 
vapour transition can be described by a single order pa- 
rameter (e.g., the density), which can take two equilib- 
rium bulk values at coexistence (the liquid density or 
the gas density). Obviously, the liquid phase will only 
wet the substrate if it is preferentially adsorbed by the 
latter. As one approaches the liquid-vapour coexistence 
from the vapour side, different scenarios are possible, de- 
pending on the substrate interactions and on the tem- 
perature: Either the liquid film covering the substrate 
remains microscopic at coexistence ("partial wetting"), 
or it grows macroscopically thick ("complete wetting"). 
The transition from partial to complete wetting can be 
first order or continuous ( "critical wetting" ) . Since crit- 
ical wetting is only expected on certain substrates at 
a specific temperature, it is rather difficult to observe 
experimentally (An experimental observation of critical 
wetting with long range forces has been reported in 0, 
and with short range forces in Q). 

Wetting phenomena are also present in alloys which 
undergo a discontinuous order-disorder transition in the 
bulk P,p|. In many cases, surfaces are neutral with re- 



spect to the symmetry of the ordered phase, but reduce 
the degree of ordering due to the reduced number of in- 
teracting neighbors. The surfaces can thus be wetted by 
a layer of disordered alloy, i.e., "surface induced disor- 
der" (SID) occurs The situation is reminiscent of 
liquid-vapour wetting; however, the underlying symme- 
try in the system restricts the possible wetting scenarios 
significantly. 

We shall illustrate this for systems with purely short 
range interactions: We consider a Landau free energy 
functional of the form 

T{m} = J dfj°dz { |(Vm) 2 + J 6 (m(r,z)) } 

d?f s (in(f,z = 0)). (1) 

Here the vector m subsumes the relevant order parame- 
ters, the z-axis is taken to be perpendicular to the sur- 
face, and dr integrates over the remaining spatial dimen- 
sions. The offset of the bulk free energy density /&(m) 
is chosen such that /(,(m&) = in the bulk. The sur- 
face contribution / s (m) accounts for the influence of the 
surface on the order parameter, i.e., the preferential ad- 
sorption of one phase or in the case of SID the disordering 
effect. In mean field approximation, the functional ([j]) is 
minimized by the bulk equation 
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dm, 
~d^~ 



difb(m), 
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which describes the motion of a classical particle of mass 
g in the external potential (— /(,(m)), subject to the 
boundary condition at z = 



g ^7 = -dif s (m) with 



dm 

dz 



V2.9/ 6 (m) (3) 
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FIG. 1. Cahn construction (schematic) for a second or- 
der wetting transition: (a) Critical wetting, (b) Partial and 
complete wetting, (c) Off bulk coexistence, approaching crit- 
ical wetting. Insets show the corresponding order parameter 
profiles. See text for more explanation. 

If the order parameter has just one component, this 
equation can be solved graphically by the Cahn construc- 
tion 0. This is illustrated in Fig. 1 for the case of a 
continuous wetting transition. The corresponding order 
parameter profiles are shown as insets. Complete wetting 
is encountered if f' s (m) crosses y/2gfb at the outer side 
of the minimum corresponding to the adsorbed phase. 
Partial wetting is found if the crossing point is located 
between the two minima (Fig. ^b). Critical wetting con- 
nects the two regimes, i.e., f' s (m) crosses \/2gJb right at 
the adsorbed phase minimum of /& (Fig. |l|a). Fig. ||c) 
shows a case where the system is off bulk coexistence. 

Now, let us consider the case of surface induced dis- 
order. Here, several equivalent ordered phases exist, 
and the ordered state breaks a symmetry. For neutral 
surfaces which do not discriminate between the ordered 
phases, fb and f s have the same symmetry. This implies 
that f s is extremal in the disordered phase (m = 0), i.e., 
\df s \ is zero at m = and thus crosses \J1gf\> at m = 
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FIG. 2. Cahn construction (schematic) for surface induced 
disorder in a system with a one-component order parameter m 
(a) at bulk coexistence and (b) off bulk coexistence. Dashed 
line shows surface term f' s (m) for critical wetting, dotted line 
for partial wetting. 



(Fig. ||). Comparing that with the scenario sketched 
above (Fig. |l|), we find that surface induced disorder- 
ing corresponds to either partial or critical wetting 
- the symmetry of the surface interactions excludes the 
possibility of complete wetting [jUj. The off-coexistence 
situation (Fig. ||b)) resembles that in Fig. |l]c). 

Alloys which exhibit surface induced disorder thus 
seem particularly suited to study critical wetting. Un- 
fortunately, the simplification due to the symmetry of 
the system often goes along with severe complications in 
other respect: Usually, one has to deal with a number of 
order parameters and other coupled fields, which inter- 
act in a way that may not always be transparent. If the 
surface under consideration does not have the symmetry 
of the bulk lattice with respect to the ordered phases, the 
interplay of order parameters and surface seg regation cre- 
ates effective ordering surface fields Jlj - H^ which may 
affect the critical behavior at the surface T fll^ , |ll| . In the 
case of a one component order parameter, such a field 
drives the system from critical wetting to partial wetting. 
When several order parameters are involved, this is not 
necessarily the case ]P6|-|l8| . More subtle effects can lead 
to surface order even at fully symmetric surfaces [|l9] ^lj . 

Experimentally, surface induced disorder has been in- 
vestigated at the (100) surface of CU3AU p^-p5[|. A num- 
ber of studies have provided evidence that the order pa- 
rameter right at the surface vanishes continuously as the 
bulk transition is approached |fj2|-p4| , and established the 
relation with the existence of a disordered surface layer of 
growing thickness [p5| . The related case of "interface in- 
duced disorder" has been studied among other in Cu-Pd, 
where the width of anti-phase boundaries was shown to 
diverge logarithmically as the temperature of the transi- 
tion to the disordered phase was approached from below 

The first simulation studies of surface induced disor- 
der in different systems have reproduced the continuous 
decrease of the surface order parameter at the bulk first 
order transition [^7 18 1(J, and the logarithmic growth 
of a disordered layer near the surface jt6|. A detailed 
study of surface induced disorder at the (111) surface 
of CuAu has been published recently by Schweika et al 
p8f . The critical behavior of various quantities has been 
analyzed, and critical exponents were found which agree 
well with the theory of critical wetting. Most notably, 
Schweika et al observe nonuniversal exponents, as pre- 
dicted by renormalization group theories of wetting phe- 
nomena |^,^l|. In contrast, Monte Carlo simulations of 
critical wetting in a simple Ising model have given results 
which were more consistent with mean field exponents 
p9| . This latter finding has intrigued theorists for some 
time, and a number of theories have been put forward 
to account for the unexpected lack of fluctuation effects 
p2| The nonuniversality of the exponents observed 
by Schweika et al seems to indicate that the fluctuations 
are restored in the case of SID. Alternatively, it may also 
stem from a competition of different length scales asso- 
ciated with the local order parameter and the local com- 



2 



position @||||. 

In the present work, we study surfaces of a bcc-based 
alloy close to the first order transition from the ordered 
DO3 phase to the disordered phase. Our work is thus 
closely related to that of Schweika et al. It differs in that 
the order parameter structure in the bcc case is much 
more complex than in the fcc-alloy: Whereas only one 
(three dimensional) order parameter drives the transi- 
tion considered by Schweika et al, we have to deal with 
two qualitatively different order parameters, which are 
entangled with each other in a rather intriguing way. In 
fact, we shall see that one of them behaves as expected 
from the theory of critical wetting, whereas the other ex- 
hibits different critical exponents, which do not fit into 
the current picture. 

A similar system has been investigated some time ago 
by Helbing et al The systems studied there were 

rather small, and a detailed analysis of the critical behav- 
ior was not possible. Helbing et al report evidence for the 
presence of a logarithmically growing disordered layer at 
the (lOO)-surface as phase coexistence was approached. 
In retrospect, this result seems surprising, since the (100) 
surface breaks the symmetry with respect to one of the 
order parameters, and we know nowadays that this nucle- 
ates an ordering surface field. In order to elucidate the 
influence of this ordering field in more detail, we have 
thus considered both the (110) surface, which has the 
full symmetry of the bulk lattice, and the (100) surface. 

Our paper is organized as follows: In the next section, 
we provide some theoretical background on the theory of 
wetting in systems with several order parameters. Sec- 
tion III is devoted to some general remarks on order- 
disorder transitions in bcc-alloys, and to the presentation 
of the model and the simulation method. Our results are 
discussed in section IV. We summarize and conclude in 
section [v| 



II. EFFECTIVE INTERFACE THEORY OF 
SURFACE INDUCED DISORDER 

A. General considerations 

We have already sketched one of the popular mean 
field approaches to wetting problems in the introduction. 
Since the bulk of the system is not critical, one can expect 
fluctuations to be negligeable for the most part. Only 
the fluctuations of the local position 1(f) of the interface 
between the growing surface layer and the bulk phase 
remain important ]3q] . As the interface moves into the 
bulk, capillary wave excursions of larger and larger wave- 
lengths become possible. These introduce long-range cor- 
relations parallel to the surface, characterized by a cor- 
relation length £11 which diverges at wetting. 

In light of these considerations, fluctuation analyses 
often replace the Landau free energy functional (|l|) by 
an effective interface Hamiltonian |pfl,b6|J37l 



H{l}/k B T = ji 



V, 



o(0}- 



(4) 



Here all lengths are given in units of the bulk correlation 
length in the phase adsorbed at the surface, the pa- 
rameter u> is the dimensionless inverse of the interfacial 
tension a 



(5) 



and the potential Vq(1) describes effective interactions be- 
tween the interface and the surface. The wetting transi- 
tion is thus identified with a depinning transition of the 
interface from the surface. 

In the linearized theory, the partition function of the 
Hamiltonian (0) is approximated by 



V{1) 



■ Jdr (VI) 2 /871 



[1+ drV (iy 



(6) 



It is convenient to switch from the real space r to the 
Fourier space q. The integration over short wavelength 
fluctuations with wavevector \q\ > A -1 , where A is ar- 
bitrary, is then straightforward: One separates I into 
a short wavelength part l(q) = l(q) 6(q — A -1 ) and a 
long wavelength part I — I — I, and exploits the relation 
Vo(l + x) = exp[cc d I dl] Vp (l) , to obtain the unrescaled 
coarse grained potential pq] 



Vx(l) 



with 



exp 



^Q 2 5m7) 



d 

dV 

l/A 1 

dq 7t 



n Jq>l/\ 



2lu ln(A/A), 



(7) 
(8) 



where A is a microscopic cutoff length. After rescaling 
r -> r/A, Vx{l) -> V\{1) = X^VxQX^), and noting that 
the roughness exponent £ is zero for capillary waves in 
d = 3 dimensions, this can be rewritten as 



VxQ) 



A 2 



2ne 



dh e -'' 2 /2«l. A V Q (l + h). 



(9) 



Renormalizing the potential Vo(l) thus amounts to con- 
voluting it with a Gaussian of width £^ A po|| , which is 
the width of a free interface on the length scale A paral- 
lel to the interface. In the case of a bound interface, a 
natural choice for A is £|| , the parallel correlation length 
of the interface. Since the remaining fluctuations after 
the renormalization should be small on this length scale, 
the procedure can be made self consistent by equating £11 
with its mean field value 



4*^(1) 



;=<*> 



= (£||Ar 2 = 1 a t A = C|| (10) 



where the average position of the interface (I) is the posi- 
tion of the minimum of V\(l). Note that the renormalizcd 
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free energy density per area £|| 2 is of order unity. The sin- 
gular part F s of the total interface free energy thus scales 
like F s oc ^ 2 . From the renormalized Hamiltonian (|J), 



H c M}/k B T 



4tt 2 



£||/A 



dq^-(q 2 + l)\l(q)\ 2 (11) 



we can now calculate the distribution probability to find 
the interface at a position h, 
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(12) 



and the joint probability distribution that the interface 
is found at h and hi at two points separated by r from 
each other. 



P {2 \h,h',f) = (S[h-l(0)] 5[h'-l(f)]) 
1 



2vr V.9(0) 2 - 9{r? 

(h-tif 



exp [• 



(ft + ti) 



'\2 



where 

ff (r) = <i(0) i(f)>~ = £ 



4(9(0) -g(r)) 4(5(0) + fl (r))- 



1 



igr/f || 



(13) 



(14) 



is the height-height correlation function of the interface 
and 



£=ff(0)«2wln(£y/A). 



(15) 



An analogous expression has been derived by Bedeaux 
and Weeks for a free liquid-gas interface in a gravita- 
tional field p9[ . In three dimensions, the height-height 
correlation function for r ^> A and £ii ^> A is a Bessel 
function K, 



g{r) = 2uKo(r/S\\). 



(16) 



We assume that the average order parameter profile 
(m(z)) is given by the average over mean field order pa- 
rameter profiles m| Ja] . e (z — I), centered around the local 
interface positions I, which are distributed according to 
the distribution function P(l). 



(m(z)) = dl P(l) m haie (z - I) 



(17) 



The distribution functions P(h) and P( 2 '(h,h',r) can 
then be used to calculate various characteristics of the 
profiles: 

For example, the effective width of the order parameter 
profile, W = 1/(2 d(m) / dz) is broadened by P(h) and 
diverges according to f40| , [4l|l 



W 2 



Wo 2 + - 6 



(18) 



where Wo denotes the "intrinsic width" of the mean field 
profile, W = 1/(2 dm bare /dz)| z=0 . 

Another quantity of interest is the layer-layer suscep- 
tibility, which describes the order parameter fluctuations 
at a given distance from the surface, 



Xzz = {(m(O)m(r)} z - (m)*}. 



(19) 



Since it has the dimension of a square length, one deduces 



immediately that Xzz scales like £ii 2 in the interfacial re- 



gion. For a more detailed analysis, we rewrite Xzz as 
Xzz = Jdf Jdh dh! m bsire (z - h) m bare (z - h!) 

x f dr {P {2) (h,ti ,r) - P(h) P(ti)}, (20) 



expand the joint probability P^ 2 \h,h' , r) in powers of 
A(r)=g(r)/a=K (r/t {l )/HZ ll /A), 



pW(h,h',r)=P(h) P( h '){l + ^ 



A(r) 



+ — + —]A W 2 + .••}, (21) 

and recall JdrrKo(r) — 1 and JdrrKo(r) 2 = 1/2. 
If the intrinsic width of the profile TibareM ^ s small 
compared to the intrinsic profile can be approxi- 
mated by a simple step profile in the interfacial region, 
TO bare( z ) = m bQ{z), where nib is the bulk order parame- 
ter. One then obtains 



Xzz = m b % 2 e-^ 2 ^ 

2w , 2w .a (£_W) 2 



(22) 



So far, these results are valid for infinite systems. The 
restriction to finite lateral dimension L affects the inter- 
face distribution P(h) ( fl2"l ) in two ways: It introduces 
a lower cutoff £\\/L in the integrals over q (e.g., (pT[)), 
and the mean position of the interface (the zeroth mode) 
is no longer fixed at the minimum of the renormalized 
potential, but distributed according to exp[— L 2 V^ (h)}. 
The width of the distribution function P(h) is now given 
by 



fll/A 



dq 



A' 



L- 



d 2 V< (h) 



dh 2 



h=(i) 



2wln(^) -wln(l + (^-) 2 ) +4vrw (^) 2 . 



(23) 



B. Bare and renormalized effective interface 
potential 

We shall now apply these general considerations to a 
specific potential Vq(1), designed to describe systems with 
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short range interactions and several order parameters and 
nonordering densities. Effective interface potentials for 
systems with two order parameters have been derived by 
Hauge |H and Kroll and Gompper [jlTj . Their approach 
can readily be generalized to the case of arbitrary many 
order parameters and nonordering densities. We choose 
the coordinate system in the order parameter and density 
space {m} such that m = in the phase which wets 
the surface, and that the coordinate axes rrii point in 
the directions of the principal curvatures of the bulk free 
energy function /&(m). Close to this phase, /& can then 
be approximated by the quadratic form 



/&(m) 



9ST 1 



- /t. 



(24) 



where /x is the field which drives the system from coex- 
istence, and the Xi have the dimension of a length. Wc 
number the coordinate axes i (i > 0) such that the Xi 
are arranged in descending order. The largest of these 
dominates the correlations at large distances and is thus 
the correlation length i.e., Xq = & = 1 in our units. 
The surface contribution has the form 



as f Q °°dl e-V-^ 2 / 2 ^Vo(l) « f^dl e -(<-«) 2 /2«i 
i.e., 

2£l - (I) < and fl/A* - (/) < (28) 

for all Xi. For Xi > 1/2, the first inequality enforces 
the second one. In a system with one order parameter, 
it leads to the well-known inequality u> < 1/2 po| , ^T| . 
As we shall see shortly, this condition is also sufficient 
to ensure the validity of the weak fluctuation limit in a 
system with several order parameters. Since to in our 
simulations turns out to be much smaller than 1/2, we 
shall not discuss the other regimes in the present paper. 

In the weak fluctuation limit, the renormalized poten- 
tial takes the form 



S ll i:A;<l/2 



a% e -i/x t (k r /xi +bi k^ e -2i +flL 

(29) 



The cutoff parameter A is of the order of the correlation 
length, A £f, = 1, and will be dropped hereafter. 



fs(m) 



hiirrii 



(25) 



C. Free energy scaling 



Following Hauge and Kroll/Gompper, we now assume 
that the actual profile from the adsorbed phase to the 
bulk phase is close to the profile of a free interface be- 
tween these two phases. Close to the surface region, we 
thus approximate the former by the test function 



mi(z) = Vi exp (z - l)/Xi 



(26) 



(at z <C I), where I denotes the position of the effective 
interface. Inserting this into eqn. (|l|) with ( pij) and (25), 
we obtain the effective interface potential 



w=E 



a,i e 



E^e-'^+^+ZiZ (27) 



for I 0, with a,i = hi.iVi and &y = ~ 9&ij I Xi)viVj. 
This expression is of course only valid for large I. Notably, 
it fails at I = 0, since the true potential Vq(1) must diverge 
there. We shall suppose that the leading term feoo = b in 
the second sum is positive and dominates over the more 
rapidly decaying terms, and disregard the latter in the 
following. 

At lo = (or in mean field approximation), the inter- 
face ist flat, and its position is given by the minimum 
of Vo(l). At nonzero to, the potential has to be renor- 
malized as described in the previous section. Now the 
renormalization is straightforward if the fluctuations are 
sufficiently small that the interface position (I) at wet- 
ting is well in the asymptotic tail of the potential (weak 
fluctuation limit). According to a criterion introduced 
by Brezin, Halperin, and Leibler |3p[| , this is true as long 



Now our task is to determine £|| self consistently by 
use of eqn. (|To|) , which will yield the scaling behavior of 
the singular part of the surface free energy, F s oc £^ 2 . 
Before generalizing to several order parameters, we shall 
briefly discuss the situation in a system with only one 
length scale Ao- The formal alikeness of the more general 
theory with this often discussed special case can thus be 
highlighted. Moreover, many of the results derived for 
one order parameter carry over directly to the case of 
several order parameters. 

In a system with one order parameter, the singular free 
energy has the scaling form 



8ttw h /($o), 



(30) 



where the scaling function f($o) depends on the dimen- 
sionless parameter 



$ = C fi (uJ - 1)/2 



a 



with C = v / (8ttw) w /2&. 

(31) 



Depending on the value of $o ; one can distinguish be- 
tween different regimes: 



$o » 1 : /(So) = 1/2 <?i(2 w $o 2 ) 

with gx(x) ss 1 + x - (2 + uo)x 2 + 
(complete wetting) 



|$o| « 1 : /(*o) w 1 - 



2 + LJ 
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(critical wetting, field like) 



(32) 
(33) 
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$0 < -1 : /($o 
with g2{x) 



($2/2)V(i-) flS (($»/2)- 1 /(M) 
3a; (2 + 7w)a; 2 



1 



2(1- 
(partial wetting) 



w) 8(1 -w) 5 



(34) 



The point $0 = is the critical wetting point. If one ap- 
proaches this point from the partial wetting side ao — > 0~ 
on the coexistence line fj, = 0, the parallel correlation 
length £|| diverges with the well-known nonuniversal ex- 
ponent 



(2iruj/b) 



-1/2(1-^) 



(-ao 



,-1/(1-") 



(35) 



and the distance between the average position of the in- 
terface and the surface diverges asymptotically like 



(0 



-(l + 2w)/(l + w) ln(-ao). 



(36) 



The relevant regime for most cases of surface induced 
disorder is however the critical wetting regime, where the 
critical wetting point is approached under a finite angle 
to the coexistence line in (ao,[J>) space. Here the parallel 
correlation length £|| scales like 



with 



'II 



1/2. 



(37) 



as fi approaches zero, the width of the interface diverges 
with 



W 2 



-wln(^), 



and its average position with 

(J)«-(w + l/2) ln(/x). 



(38) 



(39) 



These results can be used to derive the layer-bulk sus- 
ceptibility of the order parameter in the interfacial region 



Xo,(l) 



d(m ) 



Oft 



d(m ) d(l) 



dz 



1 



(0 d^i fiy/ln(fj,)' 



(40) 



In the step approximation w ]3 are (z) = mo,f>#(z), the 
layer-bulk susceptibility in the interfacial region can be 
calculated in more detail: 



X0,2 



m h 



21n/x 



(41) 



It has a slightly asymmetric peak of width £j_ at z = (I), 
the height of which scales like l//i. 

The layer-layer susceptibility could already be derived 
in the previous section. It also has a peak at the interface, 
which is however a factor of V2 narrower. Its height 
scales like 



°c £,\\ 2 /£,± 2 oc -l/(p ln(//)). 



(42) 



Next we determine the critical behavior of the order 
parameter at the surface, mo,i, 



mo i oc 



dF s 
dh i 



da 



OC jJL 



00,1 



(43) 



It will prove useful to rederive the exponent /3q i in an al- 
ternative way: The surface order parameter in mean field 
theory is given by ^bare (^) = m,bexp(— //Ao) Averaging 
the profile according to eqn. ( |T7j ) yields 

m Ql = m b (e~ l/x °) P{l) = m b e" w/Ao+ ^ /2A o. (44) 



After inserting Ao = 1 and using eqns. ( |39j ) and (|3^ 
one recovers the power law of eqn. (|4^) with the same 
exponent /?o,i- The approach has the advantage that 
it allows for a straightforward calculation of finite size 
effects on surface critical behavior: We simply replace 
the expression (|3^) for £j_ in the infinite system by eqn. 
(p3|) to obtain 



m ,i oc m b /A 1 M (8nuj ^L 1/v ^) 
with the scaling function 

M (x) = ( — 



(45) 



(46) 



We are now ready to generalize these results to the 
case of several order parameters and nonordering den- 
sities. Formally, the theory turns out to remain very 
similar. The self consistent determination of £|| leads to 
a generalized version of the scaling form for the singular 
part of the surface free energy (B0h , 



F s cx £|f 2 = 8tt W ^ /({$,}), 
where the scaling variables are 

$< = C t / ,(l-2A ! )(l-/A,)/2A, 



(47) 



(48) 



with 



As in the one-order parameter case, we have to distin- 
guish between different regimes depending on the values 
of the scaling variables. 



D. Symmetry preserving and symmetry breaking 
surfaces 



Let us first assume that the effect of nonordering den- 
sities can be disregarded {e.g., because the associated 
length scales are small, A^ < 1/2), and consider the case 
of a symmetry preserving surface. No ordering surface 
fields are then present, i.e., a,i oc hi = for all contri- 
butions i. The system is thus in a "multicritical wetting 
regime", where |$j| <C 1 for all i, and the scaling function 
can be expanded as 



G 



(49) 



am) 



J2*i$j Xj/Xi Kj(\i 



The effective interface position (/), and the correlation 
length £ii are given by eqns. (^9| ) and ((3^) as in the case 
of normal critical wetting. Hence all the results related 
to interfacial properties, such as the interfacial width, the 
interfacial layer susceptibilities etc., remain unchanged. 
In particular, the criterion for the validity of the weak 
fluctuation limit is still uj < 1/2 (from eqns. (p8[), ( |3g| ) 
and (|3£])). The surface order parameters obey the power 
law 



0£ 



rrii i oc 



oc /I 



ft 



A ' 1_ 2A; + 2^ (2Ai 



1). 
(50) 



Following the lines of (|4J) , one also obtains the finite size 
scaling function 



X +1' 



(51) 



A whole sequence of surface exponents is thus pre- 
dicted, one for each order parameter. In practice, how- 
ever, one hardly ever measures only one "pure" order 
parameter roj. Instead, one expects to observe some 
combination of contributions with different exponents 
/?t,i) which will be dominated by the leading exponent 
/3o,i = (u> + l)/2 in the asymptotic limit /i — > 0. 

The situation changes when at least one of the <Zi be- 
comes nonzero at coexistence. This is the case, e.g., at a 
symmetry breaking surface, where one or several surface 
fields become nonzero, or even at a symmetry preserving 
surface if the length scale associated with a nonordering 
density exceeds half the bulk correlation length, Aj > 1/2. 

Let a,je~ l l Xj be the leading nonvanishing term in the 
potential (f27j). As one approaches coexistence, fj. — ► 0, 
the scaling variable $j increases and one eventually en- 
ters a regime |$j| ^> 1. For negative aj, (<&j <C —1), the 
wetting becomes partial, i.e., no surface induced disor- 
dering takes place. For positive aj, (<E>j ^> 1), different 
scenarios are possible, depending on the sign and the am- 
plitude of the higher order terms cii, (i > J) in eqn. (29). 
If they are positive or sufficiently small, such that 



\a ia ~ Xj/X '\ «1, 



(52) 



the disordered phase wets the surface. The effective in- 
terface position (I) diverges asymptotically like 

(1) ^-Xj(1+uj/2X 2 ) ln(M), (53) 

the parallel correlation length scales like 



£|| = y/Xj/^uy., (54) 
and the scaling function in eqn. (j47|) takes the form 



(55) 



with 



Kj(Xi) = ^ - 1) pXjrWWA 

Xi Ai 

According to eqn. (p8|), the weak fluctuation regime here 
is bounded by u> < 2X 2 , thus encompassing the regime 
u<l/2. 

The criterion ([32]) is motivated as follows: If one of 
the higher order ai is negative and large, the interface 
potential Vf,, (I) may exhibit a second minimum closer to 
the surface, which competes with the minimum at large 
I and may prevent the formation of an asymptotically 
diverging wetting layer. The inspection of the free energy 
scaling function ([55]) reveals that the transition to such a 
partial wetting regime is appropriately described in terms 
of the combined scaling variable 



i,J 



$,:$7 Aj/Al OC a . a -W /i (A,A < -l)(l- W /(2A j A,)) 



This quantity has to be large at the point (1q where the 
one minimum of (I) splits up in two. The condition 

( f52| ) ensures that is small for all fi. 

The wetting is critical with respect to all order param- 
eters rrii with length scales Ai larger than Xj. As coexis- 
tence is approached, they vanish at the surface according 
to the power law 



9^ 



ma oc 



oc \v 



At LO , Xi 



(56) 



The finite size scaling function Mi(x) is again given by 
(|5l|), with the scaling variable x = Attuj^iL 2 / Xj. Note 
that the exponents /^i are nonuniversal even in the 
mean field limit (oj = 0). This remarkable effect has 
first been discovered by Hauge §5j| and later studied by 
Kroll/Gompper in an fee Ising antiferromagnet using a 
mean field approximation |]l7f , Monte Carlo simulations, 
and a linear renormalization group study similar to the 
one presented here fl8|] . However, (I) in this work is 
taken from eqn. ([36]) rather than determined self consis- 
tently, hence the resulting critical exponents differ some- 
what from those calculated here. As in the case of the 
symmetry preserving surface, a whole set of exponents is 
predicted by eqn. (poj). In the asymptotic limit — > 0, 
however, the surface behavior is expected to be governed 
by the leading exponent 



A, 



Aj 



2 [ Xj 



1). 



(57) 
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We have reinserted the bulk correlation length ^ = 1 
here. 

Finally, we discuss the critical behavior of the sur- 
face susceptibilities. The corresponding critical expo- 
nents can be shown to obey simple scaling laws. In the 
case of the surface-bulk susceptibility, the relation follows 
trivially: 



X»,i <* 



dfi 



OC [I 



7i,l = 1 - 



(58) 



In the case of the surface-surface susceptibility, it de- 
pends on the regime under consideration. In the "critical 
wetting regimes" discussed here, the free energy scaling 
function / can be expressed as a Taylor series in powers 
of the scaling variables $j or respectively, and one 
obtains 



dm,-,! d 2 f 
ohii oaf 



oc \x 



-7«,il 



7t,n = 1 _ 2 A,i- 
(59) 



The dominating exponents in the asymptotic limit are 
70,1 and 70,11. 



III. MODELING ORDER-DISORDER 
TRANSITION IN BCC- ALLOYS 

Fig. H shows some typical structures of binary (AB) 
bcc- alloys (e.g., FeAl Jt2|). It is useful to divide the bcc 
lattice into four fcc-sublattices a-d as indicated in the 
Figure. The phase transitions are then conveniently de- 
scribed in terms of a set of order parameters |M3j 



■01 = ( c a + Q, ~ c c - c d ) 

1p2 = (c a ~Cb + C c - C d ) 
1p3 = (Ca -C b -C c + C d ), 



(60) 



where c a denotes the composition on the sublattice a, 
i.e., the average concentration of one component A there. 
In the disordered phase, all sublattice compositions are 
equal and these order parameters vanish. The B2 phase 
is characterized by nonzero ip\, and the DO3 phase in 
addition by nonzero ip% = =b03. By symmetry, physi- 
cal quantities have to be invariant under sublattice ex- 
changes (a <-* b), (c <-> d), and (a, b) <-> (c,d). The 
leading terms in a Landau expansion of the free energy 
F thus read 



F = F + A x i\}\ + A 2 {^\ + V^ 3 2 ) + B V1V2V3 



(61) 



+c^{ + c 2 ^\ + ^ 3 4 ) + C3V2V3 2 + c^l^l + v 3 2 ), 

We point out in particular the cubic term Bipiifoipa. It 
can be read in two ways. On the one hand, it describes 
how the B2-order influences the DO3 order: The order 
parameter ipi breaks the symmetry with respect to indi- 



( a ) A2 



(b) B2 




(c) DO 3 




FIG. 3. bcc lattice with (a) disordered A2 structure (b) B2 
structure, and (c) DO3 structure. Also shown is the assign- 



vidual sign reversal of tp2 or ^3 and orients (1P2, ^3) such 
that ip2 = —sigsx(Bipi) ■03- Conversely, one can interpret 
the product 02 03 as an effective ordering field acting on 
0i. We shall come back to this point later. 

At the presence of surfaces, the situation is even more 
complicated. First, we can always expect that one com- 
ponent enriches at the surface, since there are no sym- 
metry arguments to prevent that. Even if no explicit 
surface field coupling to the total concentration c is ap- 
plied, the component which is in excess with respect to 
the ideal stoechiometry of the bulk phase ((3:1) in the 
DO3 phase) will segregate to the surface. Second, we 
have already mentioned that the Landau expansion of 
the surface free energy f s depends on the orientation of 
the surface |l^,Q. The (110) surface has the same sym- 
metry with respect to sublattice exchanges as the bulk, 
hence the Landau expansion of the surface free energy 
must have the form (pl|). In case the order is sufficiently 
suppressed at the surface, one can thus hope to find clas- 
sical surface induced disordering here. In the case of 
the (100) surface, the symmetry with respect to the ex- 
change (a, b) «-» (c, d) is broken. The surface enrichment 
of one component then induces an effective ordering sur- 
face field, which couples to the order parameter 0i [ p^[ . 
Other ordering fields coupling to 1P2 and 03 are still for- 
bidden by symmetry. The full spectrum of possible or- 
dering surface fields is allowed in the case of the (111) 
surface. 

In order to model these phase transitions, we consider 
an Ising model of spins S{ = ±1 on the bcc-lattice with 
antiferromagnetic interactions between up to next near- 
est neighbors, 

H = V ^ SiSj + aV SiSj -H^Si. (62) 

where the sum (ij) runs over nearest neighbour and ((ij)) 
over next nearest neighbour pairs. Spins S = +1 repre- 
sent A-atoms and S = — 1 B-atoms, hence the concen- 
tration c of A is related to the average spin (S) via 



c=((S) + l)/2, 



(63) 



and the field H represents a chemical potential. The pa- 
rameter a — 0.457 is chosen such that the highest tem- 
perature which can still support a B2 phase is about twice 
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as high as the highest temperature of the DO3 phase, like 
in the experimental phase diagram of FeAl. The phase 
diagram of our model is shown in Figure [|. 

The surface simulations were performed in a L x L x D 
geometry with periodic boundary conditions in the L di- 
rection and free boundary conditions in the D direction, 
varying D from 100 to 200 and L from 20 to 100. In 
order to handle systems of that size efficiently, we have 
developed Q] a multispin code f4q |, which allowed to 
store the configurations bitwise instead of bytewise jl6) . 
Our Monte Carlo runs had total lengths of up to 2 • 10 6 
Monte Carlo sweeps. 



4 




H/V 

FIG. 4. Phase diagram of our model in the T — H plane. 
Solid lines mark first order phase transitions, dashed lines 
second order phase transitions. Arrows indicate the positions 
of a critical end point (cep) and a tricritical point (tcp). 



IV. SIMULATION RESULTS 

We have studied (110) and (100) oriented surfaces at 
T = IV/ ks close to the first order bulk transition be- 
tween the ordered DO3 phase and the disordered A2 
phase. The exact bulk transition point was determined 
previously from bulk simulations by thermodynamic in- 
tegration (H, H /V = 10.00771(1) @. In the presence 
of such a high bulk field, the very top layer of a free (110) 
or (100) surface is completely filled with A particles, i.e., 
Ising spins 5 = 1. Consequently, the order parameters 
ip a and the layer susceptibilities vanish there. In the fol- 
lowing, we shall generally disregard this top layer and 
analyze the profiles starting from the second layer. 

A. (110) Surfaces - DO3 order 

We begin with a detailed discussion of surface induced 
disordering at (110) surfaces, i.e., surfaces with the full 



symmetry of the bulk. Figure |5j shows profiles of the 
order parameter of DO3 ordering per site 

One clearly sees how a disordered layer forms and grows 
in thickness as the bulk transition is approached. In order 
to extract an interface position (I) and an effective inter- 
facial width W, we have fitted the profiles to a shifted 
tanh function 

^23(n)=^ 2 b 3 lk (l + exp[-2(2-(0)/IF] )~\ (65) 

The results are shown in Figs. |] and [?]. Sufficiently close 
to the bulk transition, at (Ho — H)/V < 0.005, the data 
are consistent with the logarithmic divergence predicted 
by eqns. ([59|) and (|l8|). Intuitively, one would expect that 
an effective interface theory is only applicable if Z > W, 




10 20 30 40 

Layer 



FIG. 5. Profiles of 7/123 near a (110) surface at temperature 
T = 1 V/kB for different fields H in units of V as indicated. 
The bulk transition is at H /V = 10.00771(1). Zeroth (top) 
layer is not shown (^23(0) = 0, see text). 
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FIG. 6. Position of the interface as estimated from the fit 
(65) in units of (110) layers vs. (H - Ho)/V. 
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FIG. 7. Squared interfacial width as estimated from the fit 
(65) in units of (110) layers vs. (H — Ho)/V. Long dashed 
line shows squared interface position (Z) 2 for comparison. 

i.e., the width of the interface is smaller than the distance 
of the interface from the surface. Indeed, Fig. fj] shows 
that the logarithmic behavior sets in approximately at 
the value of H where I begins to exceed W. The pref- 
actors of the logarithms in Figs || and |7] are predicted 
to be (r/2 + u)/r)V2£b/aa in the case of (/) (Fig. |), 
and mj^/a,Q in the case of W 2 (Fig. 0), where f& is 
the bulk correlation length, ao the lattice constant, a 
factor \J~2 or 2 accounts for the distance of (110) lay- 
ers from each other in units of ao, and the parameter 
r = max(l, 2Aj/£&) depends on the length scale Xj of 
com position fluctuations (see the discussion in section 
IIP ). We shall see below that the surface data suggest 
/?l = r/2 + w(l/r - 1/2) = 0.618. Inserting this result, 
one derives 4.5 [7] < £b/a < 5.4[8] from Fig. ||, and 
£b/a > 7.8[8] from Fig. [7| These values do not agree 
with each other within in the statistical error; the in- 
terfacial width seems to decrease too fast as one moves 
away from coexistence. Yet the difference seems still ac- 
ceptable, especially considering how small the region of 
apparent logarithmic behavior is. It has been observed 
in other systems f4Sj| , that the vicinity of surfaces also 
affects the intrinsic width Wo of an interface. Moreover, 
many non-diverging terms have been neglected in eqns. 
( |39| ) and ( |l8| ) which lead to systematic errors if one is not 
close enough to H . We note that seems rather large 
for a system which is not critical in the bulk. On the 
other hand, Fig. || shows that the bulk order parameter 
ip23 decreases considerably as one approaches the phase 
transition point. This observations suggests that a crit- 
ical point is at least nearby, although preempted by the 
first order transition from the DO3 phase to the disor- 
dered phase. 

Next we consider the profiles of the layer susceptibili- 
ties of the order parameter tp23- They can be determined 
from the simulation data by use of the fluctuation rela- 
tions 



k B T 



(66) 



(67) 



where ip is the order parameter under consideration, 
N la , yet denotes the number of sites in a layer, and iV total 
the total number of sites. Fig. [s] shows that both the 
layer-bulk susceptibility Xz and the layer-layer suszepti- 
bility \zz exhibit the expected peak in the vicinity of the 
interface (eqns. ( |4l| ) and (|22])). The centers of the peaks 
can be fitted nicely by Gaussians of width £j_ and £j_ /^/2 1 
respectively, where £j_ is calculated from th e wid th W of 
the order parameter profile using £j_ = y/2/nW. The 
wings of the peaks are not Gaussian any more, but asym- 
metric - the layer susceptibilities are enhanced at the 
bulk side of the interface, and suppressed at the surface 
side. Such an asymmetry has been predicted qualita- 
tively for Xz in eqn. (fll|), but not for Xzz (cf. eqn. (p2|)). 
Even in the case of Xz, the observed asymmetry is so 
strong that it cannot be brought into quantitative agree- 
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FIG. 8. Profiles of the layer-bulk susceptibility \ z (a) 
and the layer-layer susceptibility \ zz (b) per site of the order 
parameter ^23 in units of feT, for different fields H in units 
of V as indicated. Solid line shows the fit of a Gaussian of 
width (a) £ x = (2/7r) 1/2 W / and (b) £ ± /2 1/2 to the profile 
corresponding to H — 10.004. Zeroth (top) layer is not shown 
(x(0) = 0, see text). 
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FIG. 9. Maximum of the layer-bulk susceptibility \z per 
site of the order parameter ^23 in units of fcsT vs. (Ho — H) /V 
for different system sizes as indicated. Solid line shows a fit 
to a (Ho — H)~ x behavior, and dashed line the same with 
logarithmic correction (see text). 
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FIG. 10. Maxima of the layer-layer susceptibility \ zz 
per site of the order parameter 1/123 in units of ksT vs. 
(Ho — H)/V, for different system sizes as indicated. Inset 
shows bare data, with a fit to a power law behavior with un- 
known exponent (dotted line). In the main plot, the bulk 
contribution to Xzz has been subtracted. Solid line indicates 
the slope of (Ho — H) -1 , and dashed line the whole theoretical 
prediction including the logarithmic correction. 

merit with the theory. We recall that the linear theory ap- 
proximates the capillary waves of the interface by those of 
a free interface with some suitable long-wavelength cut- 
off, i.e., they are taken to be distributed symmetrically 
about the mean interface position. The failure of the the- 
ory to describe the details of the profiles of Xz and \ zz 
presumably reflects the fact that the capillary waves are 
in fact asymmetric. Nevertheless, the main features of 
the profiles are captured by the theory. 

The centers of the peaks are slightly more distant from 
the surface than (I) in Fig. ||, but the difference is not 



significant (up to three layers at (H — H)/V = 0.0007). 
According to the theoretical prediction (|l]) and ( |42| ) , the 
heights of the peaks should diverge with 1/(H — H) with 
different logarithmic corrections. Our data are shown in 
Figs. |^ and [l0| The maxima of the layer-bulk suscepti- 
bility are best fitted by the simple l/(H — Hq) behavior, 
which the theory predicts as long as the interfacial width 
is dominated by the intrinsic width Wq. In the regime 
(Ho — H)/V < 0.005, where the capillary wave broad- 
ening of the interface becomes significant, the data are 
also consistent with th e logarithmica lly corrected version 
XT X « 1/(^0 - H)^\\n(H -H)\ (see Fig. §. 

The analysis of the layer-layer susceptibility is more 
subtle. From a double logarithmic plot of the raw data, 
one is tempted to conclude that the predicted \/(Hq — H) 
behavior is not valid; the data rather suggest a diver- 
gence with a critical exponent 0.63 (Fig. [UJ inset). How- 
ever, since we are not aware of any theoretical explana- 
tion which could motivate such an exponent, we believe 
that the apparent power law behavior over roughly two 
decades of (H — H) is most likely accidental. Look- 
ing at the values of Xzz close to the center of the slab 
(Fig. ||b)), one recognizes that the contribution of bulk 
fluctuations to Xzz is significant even close to Hq. The 
situation is complicated by the fact that the bulk fluc- 
tuations increase considerably in the vicinity of Hq, al- 
though their amplitude does not diverge. Within the 
crude approximation that the capillary waves of the in- 
terface and the bulk fluctuations are uncorrelated, one 
can subtract the latter as "background" . The thereby 
corrected data agree reasonably well with the theory, es- 
pecially when taking into account the logarithmic correc- 



tion xl 



l/(H -H)\MH -H)\ (Fig. h( 



We proceed to study the properties of the system di- 
rectly at the surface. Figure [ll] shows the order parame- 
ter V'23,1 in the first layer (recalling that the top (zeroth) 
layer is discarded) as a function of (Hq — H) for various 
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FIG. 11. Order parameter ^23,1 at the surface (first layer) 
vs. (Ho — H)/V for different system sizes L x L x D as 
indicated. Solid line indicates power law with the exponent 
/3i = 0.618. 
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FIG. 12. Finite-size scaled plot of the surface order pa- 
rameter ^23,1 vs. (Ho — H)/V for system sizes L x L x D 
as indicated. Data were scaled with exponents vn =1/2 and 
Pi = 0.618. Dashed line shows the finite size scaling function 
predicted by eqn. (51). 

system sizes. One notices finite size effects if the dimen- 
sion L parallel to the interface is small. As long as L 
is large enough, the data exhibit a power law behavior 
with the exponent pi = 0.6 18 [4]. We emphasize that 
Pi clearly differs from 1/2 here. It is close to the value 
Pi = 0.64 found by Schweika et al in their simulations of 
surface in duce d disorder in fcc-alloys ps[ . As discussed 
in section [I D , several factors may lead to such a nonuni- 
versal exponent - capillary wave fluctuations, and/or the 
presence of a length scale Aj > £&/2, which competes 
with the correlation length and would have to be asso- 
ciated with the nonordering composition fluctuations in 
the case of the symmetry preserving (110) surface. Using 
eqn. (pTj), we can derive upper bounds for the capillary 
parameter, w < 0.236, and for A,/, Aj/£b < 0.618. 

After applying finite size scaling with the exponents 
Pi and ^|| =1/2 (cf. eqn. (|45|), the data collapse onto 
a single master curve. The form of the latter can be 
calculated from eqn. (flB]), 



r/2+uj/r 

(x + ly/ 2 



(68) 



with x tx (Hq — iJ)L 1 / !y n and r — max(l, 2Aj/^), where 
the two unknown proportionality constants are fit pa- 
rameters and to = 0.236 was used (the result is only very 
barely sensitive to the choice of u>). Fig. [l^ shows that 
the data agree nicely with the theoretical prediction. 

Figure [l3| shows the layer-bulk susceptibility at the 
surface for the order parameter i{>23- According to eqn. 
( p8| ) , it should diverge with the exponent 71 = 1 — Pi = 
0.382. Indeed, the fit to our data in the region (Hq — 
H)/V < 0.02 yields 71 = 0.37[5]. In the case of the 
layer-layer susceptibility, the theory ([39]) predicts 711 = 
1 — 2Pi = —0.236, i.e., \n does not diverge at the phase 
transition. In fact, it first increases as H is approached, 
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FIG. 13. Surface layer-bulk susceptibility per site of the 
order parameter tp23 vs. (Ho~ H)/V for different system sizes 
as indicated. Solid line marks a power law with the exponent 
71 = 0.37. 

but then decreases for (Ho — H)/V < 0.02 (not shown). 
The layer-layer susceptibility at the surface here behaves 
in a similar way as observed by Schweika et al in their 
studies of surface induced disorder at the (111) surface 
of an fec-based alloy [B8fl. 
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FIG. 14. Profiles of the total composition c = ((5") + l)/2 
(top) and of the order parameter ipi (bottom) for different 
fields H in units of V as indicated. Top (zeroth) layer is not 
shown (c(0) = 1, V'i(O) = (0)). Thin dashed lines with squares 
show for comparison the profiles of ip23 from Fig. 5. 
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FIG. 15. Maximum of the layer-bulk susceptibility Xz per 
site of the order parameter in units of ksT vs. (Ho — H) /V 
for different system sizes as indicated. Solid line shows a fit 
to a (Ho — H)' 1 behavior, and dashed line the same with the 
appropriate logarithmic correction. 
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FIG. 16. Maxima of the layer-layer susceptibility \ zz per 
site of the order parameter 0i minus bulk contribution in 
units of ksT vs. (Ho — H)/V, for different system sizes as 
indicated. Solid line indicates the slope of (Ho — H)~ x , and 
dashed line the whole theoretical prediction including the log- 
arithmic correction. Inset shows bare data, with a fit to a 
power law behavior (dotted line). 



B. (110) Surfaces - B2 order 

From the results discussed so far, we conclude that the 
behavior of the order parameter 023 can be understood 
nicely within the effective interface theory of critical wet- 
ting. However, we shall see that this holds only in part 
for the second order parameter, 0i . 

Fig. shows profiles of 0i for different fields H . They 



resemble those of "023, in particular the inflection point of 
the profiles is located approximately at the same distance 
from the surface. The upper part of Fig. [li] displays pro- 
files of the total concentration c of A particles (eqn. (|63|). 
They exhibit some characteristic, ii-independent oscilla- 
tions in the first four layers, and the A concentration is 
slightly enhanced in the disordered region. However, the 
overall variation is rather small. 

The layer susceptibility profiles of the order parameter 
0i are qualitatively similar to those of 023 and not shown 
here. Fig. [Hi] demonstrates that the maximum of the 
layer-bulk susceptibility evolves with the field H as the- 
oretically predicted, X T ax oc 1/(F - H) yf\ \n(H a - H)\. 
In the case of the layer-layer susceptibility, the agree- 
ment with the theoretically expected behavior (xTz ax ~~ 
X b z f k ) oc 1/(H - H)\\n(H - H)\ is not quite as con- 
vincing, but the data are still consistent with the theory 
for (H - H )/V < 0.01 (Fig. ||). Note that the bare 
values of xTz ax would again rather suggest a power law, 
X™ aa: oc (H — H)~ 53 (Fig. [ll| inset), which is however 
most likely accidental. 

Hence the behavior of the order parameter ipi in the 
vicinity of the interface is similar to that of the order 
parameter 023 and consistent with the theory of crit- 
ical wetting. The agreement however does not persist 
when looking right at the surface. Figs, [l?] and |l^ show 
how the value of ipi in the first surface layer depends 
on (Hq — H)/V. A power law behavior is found over 
one and a half decades of (Ho — H)/V, yet the exponent 
/3i(0i) = 0.801 differs from that of 23 ,i, /3i(02a) = 0.618 
(Fig. [l7]). Moreover, the data for different system sizes 
do not collapse if one performs finite size scaling with the 
exponent v\\ =1/2 (Fig. 18(a)). The collapse is signifi- 
cantly better if one assumes that the parallel correlation 
length diverges with the exponent v\\ = 0.7 ± 0.05 (Fig. 
1(b)). 
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FIG. 17. Order parameter 0i at the surface vs. 
(Ho — H) /V for different system sizes L x L x D as indicated. 
Solid line shows power law with the exponent /3i = 0.801. 
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FIG. 18. Finite-size scaled plots of the order parameter 
ipi at the surface vs. (Ho — H)/V for system sizes L x L x D 
as indicated. Exponents are /3i = 0.801, v\\ = 0.5 in (a), and 
v\\ = 0.7 in (b). 

We have no explanation for these unexpected findings. 
The discussion in section |[| has shown that several sur- 
face exponents /^i may be present in a system with sev- 
eral order parameters. Even though we have argued that 
only the smallest exponent should survive in the asymp- 
totic limit /i — * 0, the other power law contributions 
may conceivably still dominate the behavior of certain 
quantities over a wide range of fi. However, the crit- 
ical exponent should in all cases remain invariably 
vn — 1/2. Our results seem to indicate that the behavior 
of the order parameter ipi at the surface is governed by a 
length scale, which differs from that given of the interfa- 
cial fluctuations, but which nonetheless diverges as H is 
approached. Note that v\\ w 0.7 is close to the exponent 
v = 0.63 with which the bulk correlation length diverges 
at an Ising type transition in three dimensions. Likewise, 
the exponent (3\ = 0.801 found here resembles the sur- 
face critical exponent of the ordinary transition, /3i ~ 0.8 
fl9| , j5b"f . One might thus suspect that ipi in the disordered 
surface layer becomes critical at Hq. However, such a co- 
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FIG. 19. Surface layer-bulk susceptibility per site of the 
order parameter ipi vs. (Ho — H)/V for different system sizes 
as indicated. 

incidence would seem rather surprising. Furthermore, we 
have noted earlier that the combination tp2tp3 ac ts as an 
ordering field on tpi, hence tpi cannot become critical as 
long as ?/>23 is not strictly zero. 

Figure |l^ shows the layer-bulk susceptibility at the 
surface as a function of (Hq — H)/V. It decreases as 
Hq is approached, hence the scaling relation /3i + 71 = 1 
is obviously not met for the order parameter tpn . 



C. (100) Surfaces 

Finally, we turn to the discussion of (100) surfaces. 
As already mentioned earlier, (100) surfaces break the 
symmetry with respect to the order parameter ipi, an 
ordering surface field coupling to this order parameter 
is allowed and thus usually present 0,0. This field is 
often closely related to surface segregation [ p]j40[| . In our 
case, the excess component A of the DO3 segregates in 
the surface layer and induces a staggered concentration 
field in the layers underneath, which is equivalent to ipi 
ordering. 

This is demonstrated in Fig. (2Tj. The order param- 
eters and the composition c are defined based on the 
sublattice occupancies on two subsequent layers of dis- 
tance ao/2, starting from the first layer underneath the 
surface. The top layer is again disregarded, since it is 
entirely filled with A or S = 1 . The profiles of tpi clearly 
display the signature of an additional ordering tendency 
at the surface, which in fact reverses the sign of ip\ in the 
top layers. However, the effect is rather weak and does 
not influence the system significantly deeper in the bulk. 
The profiles can be analyzed like those at the (110) sur- 
face, and mean interface positions and mean interfacial 
widths can be extracted to yield figures very similar to 
Figs. ^ and [?]. The amplitudes of the logarithmic diver- 
gences can again be used to estimate the bulk correlation 
length From the mean interface position, one 
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FIG. 20. Profiles of the total concentration of A (top, 
diamonds), of the order parameters i/>i (bottom, circles) and 
ip23 (bottom, squares) at H — 10.003 (filled symbols) and at 
H — 10.007 (open symbols). Zeroth (top) layer is not shown. 
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FIG. 21. Order parameter ip23 at the surface vs. 
(Ho — H) /V for different system sizes L x L x D as indicated. 
Solid line shows power law with the exponent f3\ = 0.61. 

calculates 4.9[7] < £fc/ao < 5.8[8], and from the interfa- 
cial width, £b/ao > 7.5 [9], in agreement with the values 
obtained for the (110) surface. Likewise, the study of 
the layer susceptibilities at the interface does not offer 



new surprises. The maxima of the layer-bulk susceptibil- 
ities for both ipi and "023 grow according to a power law 
Xz oc (Hq — if) -1 . The layer-layer susceptibility in the 
interfacial region seems to grow with a different exponent 
(~ 0.6 like in the case of the (110) surface), yet after sub- 
tracting the "background" the data are also consistent 
with the theoretically expected behavior. Last, we study 
how the surface value of the order parameter ■023 evolves 
as the transition H is approached. Fig. ^l] shows that 
it vanishes according to a power law with the exponent 
(3i = 0.61 [2], which is within the error the same exponent 
as in the case of the (110) surface. As far as the surface 
behavior of ^23 is concerned, the (100) and the (110) 
surface are thus basically equivalent. The weak ordering 
tendency of ip\ has an at most slightly perturbing effect 
on the profiles of ^23- 



V. SUMMARY AND OUTLOOK 

We have presented an extensive Monte Carlo study of 
surface induced disorder in a simple spin lattice model 
for bcc-based binary alloys. Our work complements ear- 
lier Monte Carlo simulations of Schweika et al who 
have studied surface induced disorder in fcc-based alloys 
within a similar model. Like these authors, we observe 
critical wetting behavior with nonuniversal exponents. 
We have discussed our results in terms of an effective in- 
terface model designed to describe a system with several 
order parameters. In such a complex material, nonuni- 
versal exponents may result both from fluctuation effects 
and from a competition of length scales. 

Due to the complicated order parameter structure in 
our system, however, our data could not fully be ex- 
plained within a theory which traces everything back to 
the properties of a single interface between a disordered 
and an ordered phase. The theory provides a satisfactory 
picture for the behavior of the order parameter describing 
the DO3 ordering, ^23, and in general for the structure 
in the interfacial region. However, it fails to predict the 
behavior of the order parameter of B2 ordering, ipi, di- 
rectly at the surface. Our data thus indicate that the 
fluctuations of ipi at the surface require special treat- 
ment. Parry and coworkers 34 5l[] have recently sug- 
gested an approach to a theory of wetting based on an 
effective interface Hamiltonian with two "interfaces" , the 
usual one separating the phase adsorbed at the surface 
and the bulk phase, and a second one which accounts in 
an effective way for the fluctuations directly at the sur- 
face. Our problem seems to call for such an approach. 
Unfortunately, we are far from understanding even the 
constituting elements, the fluctuations of ipi at the wall. 
We seem to observe a coupling between critical wetting 
and some kind of surface critical behavior of ipi , the ori- 
gin of which is unclear. 

Hence already our simple, highly idealized model ex- 
hibits a complex and rather intriguing wetting behav- 
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ior. In real alloys, numerous additional complications 
are present which will lead to an even richer and more 
interesting phenomenology. For example, long range in- 
teractions are known to influence wetting transitions sig- 
nificantly The effect of van-der-Waals forces on wetting 
has been investigated in detail (j]] . Van-der-Waals forces 
are important in liquid-vapour systems or binary fluids, 
but presumably irrelevant in alloys. Instead, elastic in- 
teractions caused by lattice distortions presumably play 
an important role. 

Furthermore, real surfaces are never ideally smooth, 
but have steps and islands. We have seen that the orien- 
tation of the surface affects the surface ordering. In our 
study, we did not observe dramatic differences between 
the (110) surface and the (100) surface. Nevertheless, 
we expect that the influence of the surface orientation 
on the wetting behavior can be quite substantial, e.g., in 
situations with strong surface segregation, or if surface 
orientations are involved which also break the symmetry 
with respect to the DO3 order (e.g., the (111) surface). 
Likewise, we can expect that steps and islands will af- 
fect the ordering and the wetting properties of the alloy. 
It is well known in general that the wetting behavior on 
corrugated or rough surfaces differs from that on smooth 
surfaces In addition, even a few steps or islands 

on an otherwise smooth, but symmetry breaking surface 
of an alloy can have a dramatic effect on the ordering be- 
havior, since every step changes the sign of the ordering 
surface field. 
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